Two-Band Superconductivity in MgB2 
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The study of the anisotropic superconductor MgB2 using a combination of scanning tunneling 
microscopy and spectroscopy reveals two distinct energy gaps at Ai=2.3 meV and A2=7.1 meV. 
Different spectral weights of the partial superconducting density of states (PDOS) are a reflection of 
different tunneling directions in this multi-band system. Our experimental observations are consis- 
tent with the existence of two-band superconductivity in the presence of interband superconducting 
pair interaction and quasiparticle scattering. Temperature evolution of the tunneling spectra fol- 
lows the BCS scenario jl| with both gaps vanishing at the bulk Tc. Indeed, the study of tunneling 
junctions exhibiting only the small gap (c-axis tunneling) clearly and reproducibly show that this 
gap persists up to the bulk Tc. The data confirm the importance of Fermi-surface sheet dependent 
superconductivity in MgB2 proposed in the multigap model by Liu et al. Q . 



The discovery of superconductivity in MgB2 (|] at 39K 
sparked great interest in the fundamental physics and 
practical applications of this material. There has al- 
ready been rapid progress in understanding the physi- 
cal properties of this superconductor. Specific heat mea- 
surements 1^,^ show that MgB2 is an s-wave supercon- 
ductor and the presence of the isotope effect [^j^ points 
towards phonon-mediated pairing. Tunneling and pho- 
toemission spectroscopy directly measures the supercon- 
ducting energy gap and can provide further understand- 
ing of the origin of the superconductivity in this mate- 
rial. Earlier tunneling spectroscopy measurements show 
a large spread in the gap values [p|-p^ each consistent 
with the BCS form. More recent experiments, including 
STM tunneling spectroscopy |jri|, point-contact spec- 
troscopy 1 12 1^, specific heat measurements 



and 



Raman spectroscopy |14 point towards the existence of 



two distinct gaps. This scenario has been predicted the- 
oretically by Liu et al. First principle calculations 
show that the Fermi surface of MgB2 consists of 2D cylin- 
drical sheets arising from cr antibonding states of B Pxy 
orbitals, and 3D tubular networks arising from tt bonding 
and antibonding states of B pz orbitals. In this theoret- 
ical framework two different energy gaps exist, the 
smaller one being an induced gap associated with the 
3D bands and the larger one associated with the super- 
conducting 2D bands. Furthermore both superconduct- 
ing gaps should vanish at the bulk critical temperature 
Tc. Due to this highly anisotropic band structure the 
superconducting gaps should be momentum-dependent 
reflecting the strength of the electron-phonon coupling 
of the carriers in the different bands. Up to now there 
has been no direct experimental evidence of the orien- 
tation dependence of the order parameter in this mate- 
rial. Moreover, the temperature dependence of the two 
gaps would give further insights into the nature of super- 
conductivity in MgB2. Scanning tunneling spectroscopy 
is a unique technique that allows direct measure of the 



DOS near the Fermi energy with high spatial and en- 
ergy resolution. Since the tunneling current couples only 
to bands with an appreciable Fermi velocity in the tun- 
neling direction STM can distinguish between the 2D 
bands with little or no c-axis component of Fermi veloc- 
ity and 3D bands with considerable c-axis component. 
Thus, in the clean limit (/ 27r^) the direction of the 
tunneling current selects which band and which super- 
conducting gap is probed. In this paper we report direct 
evidence of orientation-dependent double-gap structure 
in the quasi-particle energy spectra as determined from 
tunneling spectroscopy. Our data are consistent with the 
theoretical predictions of multigap superconductivity in 
the clean limit. Temperature dependent tunneling shows 
that the two distinct gaps vanish simultaneously near the 
bulk Tc. 

Compact samples of MgB2 were synthesized from 
amorphous B powder (4N's purity) and high purity Mg. 
The B powder was pressed into pellets under 6 kbar pres- 
sure. These pellets were reacted with Mg vapor at 850 
C for 2 hr in a BN container under 50 bar of Ar. Dur- 
ing the reaction the pellets broke up in pieces of several 
mm in size. The typical critical temperature of these 
gold-colored pellets is 39K. 

The tunneling measurements were performed on these 
samples using different surface preparation methods. 
One procedure consisted of chemically etching the pel- 
lets for 50 seconds in bromine (Br 1% in pure ethanol), 
then rinsing in pure ethanol and drying in N2 gas. Al- 
ternative procedures consisted of cleaving or mechanical 
polishing of the sample in an inert atmosphere. After 
these treatments the samples were mounted on the STM 
stage in He exchange gas and cooled down to 4.2 K. These 
procedures all yielded equivalent tunneling results. The 
measurements were performed with a home-built STM 
described elsewhere Q . 

The typical conductance spectra recorded on different 
grains within the sample are shown in Figure |l]. The 
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typical grain size is of the order of 50 nm and the spec- 
tra are remarkably reproducible within the same grain 
both with respect to location and tunneling resistance. 
All the spectra are normalized to the conductance value 
at -20mV. They reveal a double gap structure, flat back- 
ground, and very low zero-bias conductance with very lit- 
tle broadening other than thermal smearing. Two peaks 
are present at Vi ~ 3 mV and V2 ~ 7.5 mV, symmet- 
rically for both injection and emission of electrons and 
they stay within 10% of these values when changing be- 
tween grains and/or samples. This gap size distribution 
is in line with the proposed distribution of the gap values 
over the Fermi surface by Choi et al. fl^ ]. 

Figure H(a) is the representative surface topography 
(150 nm X 150 nm) showing several different grains. The 
roughness within an individual grain is about Inm and 
the grains are separated by steps that are 1-lOnm high. 
Each of the four grains is characterized by a tunneling 
spectrum that is consistent within a particular grain. 
Within this scanning area we focus on the tunneling spec- 
tra acquired on grains 1 and 2 (Fig. |^). The spectrum 
recorded in grain 2 shows a clear two-gap structure dom- 
inated by the high-energy gap while in grain 1 only the 
smaller gap is pronounced and the high-energy peak be- 
comes a weak satellite feature. The main characteristics 
of the spectra, including the spectral weight of each peak, 
do not change with the tunneling resistance from 0.1 to 
2 Gfl. Having this remarkable spatial reproducibility of 
the tunneling conductance we focus our attention on the 
origin of the different spectra in these grains. First, we 
point out that both spectra are in agreement with the 
multi-gap scenario proposed by Liu et al. where in the 
clean limit the impurity and surface scattering do not sig- 
nificantly influence the spectra. Second, we find a strik- 
ing similarity between the spectrum in grain 1 and the 
typical spectrum we recorded earlier on c-axis oriented 
films [|l6| . In the c-axis tunneling on MgB2 epitaxial films 
the contribution from the 3D Fermi surface is expected 
to dominate the tunneling conductance. Thus we believe 
that the distinct spectra attributable to grain 1 and 2 are 
a refiection of different crystallographic grain orientation 
with respect to the tunneling direction. This evidence 
of momentum-dependent tunneling is further supported 
by the spatial and temperature-dependent analysis of the 
spectra presented below. 

To elucidate the origin of the different spectra in the 
two grains we studied the spatial evolution of the spectra 
across the grain boundary. In Figure |^(b,c) we show the 
conductance maps at the superconducting peak values of 
Vi=-|-3.3 mV and V2=-l-7.3 mV in the region indicated 
by the square spanning the two grains. In these maps the 
bright zones mark regions where the tunneling spectra 
show higher conductance at the energy corresponding to 
Vi (Fig. |(b)) or V2 (Fig. |(c)), respectively The tran- 
sition from one type of spectra to the other evolves over 
a length scale on the order of 5 nm. We associate these 



two different types of spectra with vacuum tunneling into 
grains with two distinct orientations. The smooth tran- 
sition between the spectra across the grain boundary on 
the length scale of ^ is further evidence of the sample 
being in the clean limit. 

The distinct tunneling conductance spectra that we 
observe are consistent with the two-gap BCS model |l^] 
taking into account interband impurity scattering (see, 
e.g., 10). According to this model the partial densities 
of states (PDOS) of two bands A^i(2)(a;) 



7Vi(2)(c^) =A^i(2)(0)Rc 



^1(2) 



1(2) 



are determined by two dimensionlcss functions ui(^2){^) 
that obey equations 



MiAi — Lu + iV + iTi2 
U2A.2 — LU + iV + iT2i 



U2 — Ui 
Ui - U2 



Here Ai_2 are the gap parameters of two bands, ri2 and 
r2i are the interband scattering rates, with ri2/r2i — 
A^2(0)/iVi(0). It is important to note that strong in- 
terband scattering suppresses the superconducting tran- 
sition temperature. In particular, at r2i Tc and 
Al « A2, ST, « -7rr2i/8 (see, e. g., 

We use a linear combination of the two PDOS to obtain 
the total DOS 

Ntoti^^) = aNi{uj) + l3N2iuj) 

that we compare with the experimental tunneling con- 
ductance spectra. The PDOS' spectral weight repre- 
sented by the parameters a and (3 {a + l3—l) is associated 
with the relative orientation of the crystallographic axis 
of a grain with respect to the tunneling barrier. Thus 
in the fits shown below the parameters a and (3 are kept 
constant for all spectra obtained on a particular grain at 
different temperatures. 

The linear combination of the two PDOS shown above 
accounts for all major features in our tunneling spectra 
including the temperature dependent data. The evolu- 
tion of the tunneling spectra in the temperature range 
between 4.2K and 42K are displayed in Figure ^. The 
two-gap signature in the spectra of grain 1 becomes 
hardly distinguishable with increasing the temperature 
(Fig. ^(a)). The small gap shows a rapid decrease in the 
peak height and the conductance peaks of the two gaps 
merge into a broad peak at T«10K. In order to sepa- 
rate the temperature evolution of each PDOS we used 
the temperature dependence of the spectrum on grain 2 
(Fig. |^(b)) where only the small gap is present. Each 
curve was normalized to the conductance value at -20 
mV and then compared with the theoretical curve at the 



2 



same temperature. The theoretical curves were obtained 
by fixing the interband scattering rates ri2, and the 
smearing parameter F at values optimized for the 4.2 K 
conductance spectra. In all the fits we put Ti2= 0.73 
since their ratio is proportional to the ratio of the 
density of states at the Fermi level of the 3D and the 2D 
bands ||l8|,^ . The theoretical curves reproduce very well 
the shape inside the gap but not the peak heights. By 
adjusting the parameters F12 and F we could reproduce 
very well the peak heights but the zero bias value would 
be higher than the experimental one. Fits to the two set 
of parameters have been used to get an estimate of the 
error bars for the two gap values: Ai= 2.3±0.2 mV and 
A2= 7.1±0.2. A very good fit of both peaks' height and 
zero bias value can be obtained only for interband scat- 
tering values that are very large ('^2.5 mV) and values of 
the 3D gap is extremely small. According to jl^] such a 
large interband scattering value should lower the critical 
temperature by ^--^9 K which contradicts our experimen- 
tal findings. The set of temperature dependent tunneling 
spectra clearly shows that both superconducting gaps ex- 
ist up to 37K. 

Next, we turn our attention to the spectral weights of 
the two PDOS to the total DOS. The particular choice 
of the two grains in Fig. ^ was based on the peculiarity 
of the conductance spectra of these crystallites. They 
closely represent the two extreme limits of partial contri- 
bution by each PDOS that we have observed so far on a 
large number of pellets and thin films. The theoretical 
fits obtained by our model in these two limits (Fig. ^,b) 
show that the contribution of the 3D PDOS varies from 
nearly 100% (Fig. |b and down to 87% (Fig. 

Therefore it appears that momentum averaged STM tun- 
neling conductance is dominated by the 3D band in this 
momentum dependent two-band tunneling process. This 
could be due to the specific geometry of the tunneling 
cone emanating from the STM tip and its projection on 
the highly anisotropic Fermi surface of the MgB2 crys- 
tallite. Obviously, the relative weights of the PDOS are 
temperature independent and the parameters a and /3 
remain fixed when the temperature dependent fits of the 
tunneling spectra are obtained. 

The values of Ai(T) and A2(T) are extracted from the 
theoretical curves and the results are reported in Fig- 
ure ||(c). The high-energy gap follows a BCS behavior. 
The smaller gap remains constant up to 20K and follows 
the BCS predictions with a critical temperature far be- 
yond the one expected for such a small gap (Tc «14K) 
with a ratio 2A/kTc wl.4 similar to earlier findings from 
point spectroscopy p3| . Such behavior clearly demon- 
strates that the gap in the 3D band is mainly induced by 
the interband pairing interaction rather than produced 
by pairing interaction inside this band, in agreement with 
theoretical expectations This experimental evi- 

dence directly excludes the possibility that the grains or 
their surfaces are associated with different stoichiometry. 



The high-energy gap follows a BCS behavior with a ra- 
tio 2A2(0)/fcrc ~ 4.3, significantly larger than the BCS 
value of 3.52. There are two possible sources for this in- 
crease. Usually such enhancement is attributed to strong 
coupling effect. However, this ratio also increases due to 
the presence of additional superconducting bands even in 
the weak coupling limit. In the case of Ai ^ A2 one can 
derive a simple formula from the two-band BCS theory 

2A2 / A^iA? A2\ 

Using A''2 — 0.73Ni from |^ and experimental values 
of the gaps, we estimate from this formula 2A2/TC w 
4.1. Therefore, we conclude that the main part of the 
2A/Tc increase originates from the multiband nature of 
the superconductivity in this compound. 

Finally, we would like to point out that in the case 
of two-band superconductivity in the dirty limit {I <^ ^) 
the conductance spectra are very broad and similar to 
those reported earlier These spectra are mostly ob- 
tained on pellets without any surface treatment. Im- 
purities near the surface may provide sufficiently strong 
interband scattering to smear out the two-gap features 
in this case. 

In conclusion, the existence of two gaps has been ob- 
served by tunneling spectroscopy in compact pellets of 
MgB2. Tunneling spectra show different ratios between 
the two gap heights on different grains within the sam- 
ple. This result can be interpreted in terms of different 
tunneling directions with respect to the crystallographic 
orientation of the grain, therefore supporting the two- 
gap scenario (|Jl^. Different tunneling direction should 
probe the two bands with different weights. Moreover, 
the two gaps both vanish at a temperature close to the 
bulk value. The observed temperature independence of 
the small gap up to T^v 20 K implies that no transition 
temperature corresponding to the BCS ratio 2Ai=3.52 
kTc is present. This suggests that the coupling between 
the two bands is moderately strong and that the super- 
conductivity is dominant in the 2D sheets. 

This work was supported by US DOE Basic Energy 
Science - Material Science under contract No. W-31-109- 
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FIG. 2. (a) Topographic image of the scanning area 150 
nm X 150 nm acquired in the constant current mode (sample 
bias V= -20 mV, 1= 200 pA). The conductance spectra and 
the current-voltage characteristics on two different grains are 
shown on the right. The tunneling resistance is 0.2 Gfl. (b) 
and (c) Conductance maps of the scanning area 30 nm x 30 nm 
outlined in image (a) by the square across the grain boundary. 
The conductance maps have been recorded at Vi=3.3 mV (b) 
and V2=7.3 mV (c) respectively. 
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FIG. 1. Tunneling conductance spectra recorded on differ- 
ent grains at T=4.2K. The tunneling resistance is 0.2 GQ. 
The spectra have been normalized to the conductance value 
at high voltage. 
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FIG. 3. Temperature evolution of the two tunneling spec- 
tra together with the theoretical curves. For all the curves 
ri2=0.1 mV, r2i=0.14mV, r=0.11mV (values optimized to 
the experimental curves at 4.2 K). In (a) the curves are re- 
produced using a= 0.87, /3=0.13. In (b) a= 1, /3=0. In (c) 
the gap values axe extracted from the theoretical curves and 
are plotted as a function of the temperature together with the 
BCS A(T). 
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